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Abstract 



We re-analyze proton beam dump data taken at the U70 accelerator at IHEP Serpukhov 
with the z^-calorimeter I experiment in 1989 to set mass-coupling limits for dark gauge 
forces. The corresponding data have been used for axion and light Higgs particle searches 
in Refs. [1,2] before. We determine new mass and coupling exclusion bounds for dark gauge 
bosons. 



1 Introduction 



Long range forces based on a U(l) gauge symmetry beyond those of the SU (3) c x SU (2) l xU(1)y 
Standard Model may exist yet unnoticed if their coupling to ordinary matter is very weak [3]. 
Symmetries of this kind are discussed in various extensions of the Standard Model, see the 
surveys [4-6] . A new U{ 1 ) gauge boson 7' with masses my in the Me V-Ge V range extends the 
Lagrangian of the Standard Model £sm to [4, 7] 

C = C SM - -X^ + ^X^F 1 " + e*£h^X» + ^X^X" . (1) 

Here X^ denotes the new vector potential and X^ v = d fl X u — d u X^ the corresponding field 
strength tensor, with F^ u the U(l)y field strength tensor. The mixing of the new U(l) and 
U(l)y of the Standard Model is induced by loops of heavy particles coupling to both fields [3,6]. 
The field X^ is assumed to couple minimally to all charged Standard Model fermions i/j, with 
effective charge e^e, where is the fermionic charge under t/(l)y. For the generation of the 
mass term we assume the Stueckelberg formalism [8] , as an example. Possible other mechanisms 
consist in technicolor or spontaneous symmetry breaking. The latter ones would lead to more 
terms in (1). The parameter e denotes the mixing parameter of the two U(l) groups and may 
take values in the range e ~ 10~ 23 — 10~ 2 , depending on the respective model, cf. [4]. 

Dark U(l) gauge forces contribute to the anomalous magnetic moments of the electron and 
muon. Potential signals may be measured from T(3S) decays. The 7' particles may be created 
in electron- and proton beam dumps. So far signals of these particles have not been detected 
leading to various exclusion bounds in the my — e plane in the range of e G [5 x 10~ 9 , 10 -2 ] and 
a series of mass regions in my e [2 m e ,3 GeV], cf. e.g. Refs. [4-6] and [9]. 

In the present note we derive new exclusion bounds on dark 7' bosons using proton beam 
dump data at p ~ 70 GeV. These data have been used in the in axion [10] and light Higgs boson 
searches, cf. [1,2,11] in the past. We first describe the production process and the experimental 
facility, and then derive new mass and coupling bounds. 



2 The Production Process 

The abundant production of 7r° mesons in proton beam dumps leads to a production rate of 7' 

/ m 2 \ 3 

a(pp j'X) = 2e 2 1 f) Br(7T° 77) <r(pp -> tt°X) , (2) 

through 7T° — > 77, with the mixture 7 = (7 + ej')/(l + e) and Br (tt° 77) = 0.98823 ± 
0.00034 [12], neglecting contributions in higher powers of e. Here the phase space factor of the 
2-particle decay has been accounted for [13]. The mass range of the produced 7' is limited by 
my < m n o = 134.976 MeV [12]. In principle this range may be expanded to higher meson 
decay thresholds including r],p ± ,u, and rj' production, since the decay spectra of these particles 
contain a large fraction of photons, see also [7, 14] . However, one has to know the corresponding 
differential meson production spectra for pp resp. pA scattering in detail, which to our knowledge 
have not been measured in the energy region under consideration. 

We consider the collision process pFe — > ir°X at a momentum of the incoming proton of 
p = 68.6 GeV. The differential scattering cross sections for the reactions pp — > ir ± X for values 
of p± < 1 GeV were measured [15]. One may use the representation 



E 



d 2 a(pp -> 7T°X) _ I 
dxpdp 2 , 2 



^d 2 a(pp ->• tt + X) ^d 2 a(pp -> ir~X) 
dxpdp\ dxpdp\ 



(3) 
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for the invariant cross sections. The differential cross sections were parameterized by 

E d2aipP ^f X) = A± exp(B ± \x F \ + C±4) exp(D ±P± + E ± p\) (4) 

in [15]. Here ii? = Pl/Pl*** denotes Feynman-x and p± the transverse momentum. The normal- 
izations obey A+/A- = 2.16 ± 0.24 and the other parameters are 





+ 




B 


-5.21 ± 1.03 


-9.51 ±0.21 


C 


-1.80 ±2.62 


+2.14 ± 1.74 


D 


-1.80 ±0.31 


-1.22 ±0.38 


E 


-4.26 ±0.48 


-4.44 ± 0.58 



Here we refer for the p± distribution to the values of (p±) = 0.2 GeV. The resulting distributions 
in p\ and xp are shown in Figures 1,2. 




The actual beam dump experiment used an iron target. The inclusive iron-proton cross section 
is related to the proton-proton cross section from (2) by 

a(pFe -»■ 7T°X) = A a{XF) a(pp tt°X) , (5) 

with A = 56. The A a ^ XF ^ dependence of the 7r°-pro duct ion cross section on nuclei as of iron can 
be parameterized by 

(a{x F )) ~ 0.55 , (6) 
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based on the compilation given in [16]. The inclusive cross section at 69 GeV was measured 
with [17] 

a(pp -»■ 7T°X) = 74 ± 1.0 mb . (7) 




Figure 2: The p 2 ^ and xf distribution of the 7r° production spectrum in pFe collisions at p = 68.6 GeV. 

We generate the 7' particle production in the pp center-of-momentum system and boost to 
the laboratory system then. We investigate the mass range of 7' particles above the 2-electron 
threshold of 1.022 MeV up to m n o. A cut {px/PLjiab < tan(9 max ) is applied to ensure that the 
produced particle reaches the fiducial volume of the detector (see below). The effect of this cut 
is illustrated in Figure 3 which shows the fraction of particles which pass the fiducial volume cut 
with respect to all produced particles per energy bin. Results for 7r° are given as well as for 7' 
of various masses. The kinematic factor (1 — my/m^) 3 from (2) is included for the 7'. 



> 




(GeV) 

Figure 3: Fraction of produced 7' particles which reach the detector as function of their energy in the 
laboratory frame. Colors: 7' with masses between and 120 MeV in steps of 20 MeV (from top to 
bottom). The black line shows the corresponding 7T° for comparison. 
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The only relevant decay channel in the considered mass range is 7' — > e + e . The inverse live 
time of the 7' particle is then given by [7] 

. 1 1 „ I W" / 2ml\ 

r(V)- 1 = rQE om-,,^l - ^ [l + ^ j . (8) 

Potential electron-pairs from 7' decay manifest as electromagnetic showers in the detector used. 
The mass and coupling limits on the 7' bosons given below are derived from the observed rate 
of these showers over the background, cf. [1]. 



3 The Experimental Setup and Data Taking 

The beam dump experiment was carried out at the U70 accelerator at IHEP Serpukhov during a 
three months exposure in 1989. Data have been taken with the I/-CAL I experiment, a neutrino 
detector. All technical details of this experiment have been described in [1] and a detailed 
description of the detector was given in [18]. Here we only summarize the key numbers which 
are crucial for the present analysis. 

The target part of the detector is used as a fiducial volume to detect the decays of the 7'. It 
has a modular structure and consists of 36 identical modules along the beam direction. Each of 
the modules is composed of a 5 cm thick aluminum plate and a pair of drift chambers to allow 
for three dimensional tracking of charged particles. 

For the beam dump experiment a fiducial volume of 30 modules with a total length of I = 23 m 
is chosen, starting with the fourth module at a distance of /dump = 64 m down-stream of the 
beam dump. The lateral extension of the fiducial volume is 2.6 x 2.6 m 2 . In the following we use 
conservatively a slightly smaller fiducial volume, defined as a cone pointing to the beam dump 
with a ground circle of 2.6 m in diameter at the end of the fiducial volume, i.e. at a distance of 
87 m from the dump. This leads to the following simple fiducial volume cut: 

(p±/p L )iab< 1.3/87 = 0.015 . (9) 

For particles which traverse the fiducial volume, the decay probability Wdec for decays 7' — > e + e~ 
is then given by 



Wdec = exp 



/dump TTly' 

cr(7') p 



1 — exp 



I m^i 
ct(j') p 



(10) 



with c the velocity of light, my and p are the mass and momentum of the 7'. 

During the three months exposure time in 1989 1.71 x 10 18 protons on target had been 
accumulated [1]. The signature of event candidates from 7' — > e + e~ is a single electromagnetic 
shower in beam direction. This signature is identical to the one from the axion or light Higgs 
particle decay search which was performed in [1]. For energy deposits above 3 GeV in the 
detector the reconstruction code can distinguish electromagnetic from hadronic showers very 
well. Therefore the following final cuts for the selection of isolated electromagnetic showers had 
been chosen [1] : 

• A minimal electromagnetic shower energy of E c \ m > 3 GeV; 

• A maximal hadronic shower energy of -E^ad < 1.5 GeV; 

• A maximal angle with respect to the beam direction of O cmi < 0.05 rad. 
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From the total data sample of 3880 reconstructed events, 5 pass these cuts. Background estimates 
from the simulation of and v e interactions in the detector account for 3.5 events. The Poisson 
probability to observe 5 or less events for an expectation of 3.5 events is 86%. Data are therefore 
compatible with the simulated background from conventional neutrino interactions. 

4 Search for decays 7' — > e + e~ 

Signals from 7' — > e + e~ decays pass the cuts mentioned above with an energy independent 
efficiency of e = 70% if the true energy of the decaying particle is above 3 GeV. The total 
number of expected signal events can therefore be calculated as 

at at &(pFe 7'X) forward 

N sig = N tot x — x w; dec x e , (11) 

a(pp —> X) 

with iVtot the total number of protons on target during the exposure time. The index 'forward' 
indicates the application of the fiducial volume cut. The dependence of iV si g on my and e is 
shown in Figure 4. 




Figure 4: Detected 7' events. Color scale in log 10 from 10 8 events (red) to one event (dark blue). 

For my ~ 2 MeV and e ~ 2 • 10~ 4 the decay probability reaches its maximum of 11% and more 
than 10 8 signal events would be expected in the detector. For larger e the decay length decreases 
exponentially and at some point all 7' particles decay before reaching the detector. For smaller 
values of e the decay length increases. At e < 10~ 6 most of the particles pass the detector without 
decaying. For even smaller values of e the event rate of both production and decay probability 
decrease proportional to e 2 . The dependence of the iso-event number lines on my is governed 
by the boost factor p/my. The sensitivity is kinematics limited to values 

For 10.6 expected events the Poisson probability to observe 5 or less events is less than 5%. 
With a background expectation of 3.5 events we can therefore exclude a signal contribution of 
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7.1 events at 95% CL. The corresponding exclusion region is shown as the red line in comparison 
with limits from other experiments in Figure 5, see also [4,6,7]. 



CO 




10" 2 10" 1 1 

m y '[GeV] 

Figure 5: Comparison of the present exclusion bounds (red line) with other limits from the measurement 
of the anomalous magnetic moments a e and a M [19], T(3S) decay [20], the beam dump experiments 
E137, E141, E774 [21-23], and supernovae cooling [4,24]. We indicate the prospects for LSND [7,25] 
(open grey-bounded area), and the DAMA/LIBRA region (open orange bounded area) [26]. The limits 
for e > 1CP 7 have been taken from Ref. [6]. 

At large values of e studies of the anomalous magnetic moments of the muon and electron [19] 
and of rare decays of heavy mesons [20] put stringent limits. For 10~ 3 < e < 10~ 7 beam dump 
experiments [21-23] give the best sensitivity. For even smaller values of e limits can be derived 
by studying the dynamics of supernovae cooling [24]. For completeness, the prospects for the 
sensitivity of a reanalysis of LSND data [25] is also shown, as proposed in [7]. 

The present analysis is sensitive to the region 10~ 3 < e < 10~ 7 and the exclusion region 
largely overlaps with the one from E137. However, we are sensitive to larger values of e for the 
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same values of my. In particular, a new region around e ~ 10~ 5 and my « 50 MeV is explored 
here. This region is part of a band ( shown in orange) which correlates e and my within certain 
supersymmetric theories [27]. Within these models heavy dark matter candidate particles pick 
up an U(l)' charge and can scatter elastically through 7' exchange. This would allow to explain 
the annual modulation of the DAMA/LIBRA experiment [26]. For a more extensive discussion, 
see [4]. 

5 Conclusions 

We have re-analyzed proton beam dump data taken at the U70 accelerator at IHEP Serpukhov 
with the //-calorimeter I experiment in 1989 to set mass and coupling limits for dark gauge forces, 
searching for electromagnetic signatures according to the decay 7' — > e + e~ . The analysis extends 
the region excluded by former experiments in the mass region my £ [0.03 GeV, m no ] towards 
larger values in the mixing parameter e £ [2 x 10~ 6 ,2 x 10~ 5 ]. A lower part of the anticipated 
DAMA/LIBRA region is excluded. At lower values of e ~ 10 -7 a smaller region of masses in 
the range my £ [0.03,0.1] GeV is excluded beyond the bounds given by E137 [21]. In future 
experiments signals from dark gauge forces will be searched for in the yet unexplored regions 
shown in Figure 5, see e.g. Ref. [6] for proposals. 
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